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We study the dynamics of the chiral phase transition expected during the expansion of 
the quark-gluon plasma produced in a high energy hadron or heavy ion collision, using 
the 0(4) linear sigma model in the mean field approximation. Imposing boost invariant 
initial conditions at an initial proper time t and starting from an approximate equi- 
librium configuration, we investigate the possibility of formation of disoriented chiral 
condensate during the expansion. In order to create large domains of disoriented chiral 
condensates low-momentum instabilities have to last for long enough periods of time. 
Our simulations show no instabilities for an initial thermal configuration. For some of the 
out-of-equilibrium initial states studied, the fluctuation in the number of particles with 
low transverse momenta become large at late proper times. 

1. MODEL AND APPROXIMATIONS 

The possibility of producing large correlated regions of quark condensate < q~iqj > 
pointing along the wrong direction in isospin space was proposed [1] to explain rare events 
where there is a deficit or excess of neutral pions observed in cosmic ray experiments 
[2]. The idea that such disoriented chiral condensates (DCC's) can be formed has been 
the source of several experimental proposals [3] in high energy collisions. It has been 
recognized that nonequilibrium dynamics can yield regions of DCC [4-6] (see Gavin and 
Asakawa contributions in these proceedings). 

There are clearly two important questions that the theoretical models should aim to 
answer. The first one is to determine whether during the evolution that follows the 
collision there are instabilities affecting the fluctuations, in which case there is a chance 
of the correlations growing. The second question is, assuming the instability exists, can 
the correlated domains grow large enough so that many pions can be emitted from each 
domain making the detection of DCC's possible. 

To address these questions we employ the 0(4) linear sigma model, using boost invariant 
initial conditions at an initial proper time t , in the mean field approximation. This 
model seems to have the essential attributes of being simple but realistic enough: it has 
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the correct chiral symmetry properties, it also appropriately describes the low energy 
phenomenology of pions and dynamics then incorporates a cooling mechanism that may 
lead to low momentum instabilities. Moreover, the mean field approximation improves 
the classical approximation, because it takes into account quantum fluctuations that can 
not be neglected when the system cools down, and because it enables us to estimate the 
number of pions and its fluctuations at late times when the condensate settles down near 
to its vacuum expectation value. 

In this model the fields are organized in an 0(4) vector (</>i, </> 2 , </>3, $4) = (cr, 7Ti, 7r 2 , 7r 3 ) 
and we use the following set of coordinates 

t-(* 2 -* 2 ) 1/2 , i/^log(^), x ±> (1) 

Z t — Z 

where t is the proper time and r] is the spatial rapidity. The equations of motion are 
given by 

r-'dr rd T &(r) + x(r) &(r) = H6 iu (2) 
X(r) = X(-v 2 + 2 (r)+ < <p 2 ( V , x ± , t) >]), (3) 

i=l 

(r- 1 d T rd T - T - 2 d 2 v -dl + x(r))^ V ,x ± ,r) = 0, (4) 

where % = X(z^iL\ 4>i4*i — v 2 ) with N = 4. At equal proper times the quantum field <pi(x, t) 
defined here satisfies 

< <p(r]x,x ± ,T) <p(ri y ,y ± ,T) >= G c (rj x ,x ± ,r] y ,y ± ;T), (5) 

where G c (x,y) = jjJ2f =1 < (j>i(x)(j)i(y) > —<j)i(x)(j)i(y). We expand the field ip in an 
orthonormal basis 

<p(ri, x ± , t) = J dk(exp(ik ■ x)/ fc (r) a k + h.c), (6) 

where k ■ x = k v r] + k±x±, dk = dk T1 d 2 k±/(2ir) 3 , and the mode functions /^(t) evolve 
according to: 

h = -(^| + V 2 + X(r) + = -ulfk\ V = \k±\. (7) 

A dot here denotes the derivative with respect to the proper time t. The Fourier transform 
of the expectation value < (c^) 2 (x,t) > can be expressed in terms of the mode functions 
fk and of the distribution functions 

n k =< a[a k >, g k =< a k a k >, gl =< a[a[ > (8) 

which entirely characterize the initial state of the quantum field. 

In order to solve equations (2), (3) and (7) as an initial value problem, we need to fix at 
t the mean values (j)i k , <j> ik , n k , g k and g\. In the boost invariant case ^(t) = ^(t) = 
for k ^ 0. 

The complexity of the processes we study suggests that our initial configuration is not 
a pure state but a mixed ensemble that should be described by a density matrix. In the 



3 



mean field approximation the most general density matrix for each mode & is a gaussian 
with five parameters 1 associated with the above initial conditions and it can be written 
as 

p(&k, </>'ik) = J ' exp{i<f> lk (</)'* k - 4>* k ) - \{K lk - AG lk T, 2 k )\</)' lk - (j) lk \ 2 

-iH lk [\(j) lk - <j) lk \ 2 - \<j)' lk - fa k \ 2 ] - —^-\<j) lk + <j)' lk - 2fa k \ 2 }, (9) 

where 

G lk = < (j) lk (j)* k > -<f> lk (f>* k ; K ik =< cj) lk <j)* k > -(f> lk cj) lk (10) 
4E lfc G lfc = < {(j) lk , <j)* k } > -{(j) lk , <j) lk }. (11) 

Pure density matrix satisfies the condition p = K — 4GS 2 + 1/4G = (no (^{(tyft^k) 
mixed term). This subset corresponds to squeezed-coherent states [7]. For E = and 
Gi k = 1 1 \f2sjOk, the density matrix operator becomes p c = \c >< c\ where |c > is a 
coherent state. One can make a direct connection between our calculation and earlier 
work [4], which used classical approximation to the same model. The equilibrium density 
matrix of a free theory is Gaussian and in the high T limit (/3o>fc << 1) has the following 
representation: 

?=II / ^4 fc exp{-^(|^| 2 + (^ 2 +m 2 )|^| 2 )}^(^ ,y. (12) 

ik 

Therefore, solving the mean field equations with an initial thermal homogeneous density 
matrix is equivalent to solving the semi classical equations for a thermal ensemble of 
inhomogeneous "classical" initial conditions. 



2. RESULTS 

Fig. 1 shows the mean density of particles < n{jp) >= p J dk v < a k (t)a,k(t) > /(2ir) 3 = 
p J dk v n k (t) I (27r) 3 and the fluctuation An(p) = pj dk v < (n k (t)— < n k (t) >) 2 > /(2tt) 3 
at the initial proper time t = lfm/c and at t = 10, 20/m/c. The figures on the left 
correspond to an initial thermal configuration. In this case % does not become negative 
during the evolution so no instabilities develop. The density < n(p) > hardly changes 
during the evolution, as does its fluctuation. The fluctuation at large proper times and at 
the initial proper time are of the same order of magnitude. Therefore, observing a DCC 
for an initial thermal ensemble is unlikely. In the right column we generated instabilities 
by choosing an initial ensemble that is far from equilibrium for which we modified only 
the initial condition for a from to —1. In this case % becomes negative at early stages of 
the evolution, and an exponential growth of the low momenta modes occurs. The density 
< n(p) > is strongly amplified in the region of small transverse momenta by acquiring 
energy from the a field . The fluctuation of the number density also grows substantially for 

■"^We note that for a gaussian density matrix the requirement of isospin invariance leads immediately to a 
factorization form p = Y\i Pi- Due to the factorizability of the density matrix, if one species i has a large 
average < N l >, the whole density matrix will also have a large < N > and vice versa. This is also true 
for the variance < (N - N) 2 >. 
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Figure 1. The mean density of particles < n{jp) > and the fluctuation An(p) at t = 1/m/c 
and at t = 10, 20/m/c. The left column is for an initial thermal distribution where 
n\ = l/(exp/3a; fc - 1), ct(t = 1) = a eq (T = 200MeV), vf l (l) = ?P(1) = £(1) = 0. The 
right column has the same initial conditions as in the left one, but for cr(l) = —1. 



small transverse momenta, and the observation of DCC is more likely for this "peculiar" 
ensemble. For similar out-of-equilibrium ensembles we found no instabilities if we start the 
simulations at proper times greater than If raj c. This implies that if the expansion rate is 
too slow or if we start with lower energy densities even for out of equilibrium ensembles, 
DCC will hardly be produced. On the other hand populating at t the low momenta 
fluctuations in a thermal like distribution with a temperature T\ and populating the 
rest of the higher momenta in a thermal like distribution with a lower temperature T 2 , 
we found that instabilities occur even if we don't put kinetic energy in the condensates. 
Such an initial configuration is close in nature to the quench approximation. 
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